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In this talk we discuss the BFKL Pomeron Calculus and its interrelation with the 
colour dipole approach. The two key problems we consider are the probabilistic 
interpretation of the BFKL Pomeron Calculus and the possible scenario for the 
asymptotic behaviour of the scattering amplitude at high energy in QCD. 



1. Introduction 



The simplest approach that we can propose for high energy interaction is 
based 1,2 on the BFKL Pomeron 3 and reggeon-like diagram technique for 
the BFKL Pomeron interactions 4 < 5 ' 6 ' 7 . This technique which is a general- 
ization of Gribov Reggeon Calculus 8 , can be written in the elegant form 
of the functional integral (see Ref. 5 and the next section); and it is a 
challenge to solve this theory in QCD finding the high energy asymptotic 
behaviour. However, even this simple approach has not been solved dur- 
ing three decades of attempts by the high energy community. This failure 
stimulates a search for deeper understanding of physics which is behind 
the BFKL Pomeron Calculus. In particular, at the end of the Reggeon 
era it was understood 9 ' 10 > n that the Reggeon Calculus can be reduced 
to a Markov process 12 for the probability of finding a given number of 
Pomerons at fixed rapidity Y. Such an interpretation, if it would be rea- 
sonable in QCD, can be useful, since it allows us to use powerful methods 
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of statistical physics in our search for the solution. 

The goal of this talk is to consider two main problems: the probabilistic 
interpretation of the BFKL Pomeron Calculus based on the idea that colour 
dipoles are the correct degrees of freedom at high energy QCD 13 and the 
possible solution for the scattering amplitude at high energy. Colourless 
dipoles play two different roles in our approach. First, they are partons 
('wee' partons) for the BFKL Pomeron. This role is not related to the large 
N c approximation and, in principle, we can always speak about probability 
to find a definite number of dipoles instead of defining the probability to 
find a number of the BFKL Pomerons. The second role of the colour dipoles 
is that at high energies we can interpret the vertices of Pomeron merging 
and splitting in terms of the probability for two dipoles to annihilate in one 
dipole and of the probability for the decay of one dipole into two. It was 
shown in Ref. 13 that P — > 2P splitting can be described as the process 
of the dipole decay into two dipoles. However, the relation between the 
Pomeron merging (2P — ► P) and the process of annihilation of two dipoles 
into one dipole is not so obvious and it will be discussed here. 

Our presentation is based on Refs. 14 ' 15 ' 16 ' 17 ' 18 and we would like to 
thank Michael Lublinsky and Alex Prygarin for their contributions and for 
a great pleasure to work with them. 

The outline of the talk looks as follows. In the next section we will 
discuss the BFKL Pomeron Calculus in the elegant form of the functional 
integral, suggested by M. Braun about five years ago 5 . We will show 
that the intensive recent work on this subject 19 ' 20 > 16 confirms the BFKL 
Pomeron Calculus in spite of the fact that these attempts were based on 
slightly different but not more general assumptions. 

In the third section the general approach based on the generating func- 
tional will be discussed. The set of equations for the amplitude of n-dipole 
interaction with the target will be obtained and the interrelations between 
vertices of the Pomeron interactions and the microscopic dipole processes 
will be considered. 

The fourth section is devoted to the toy model which simplifies the QCD 
interaction and allows us to see the main properties of high energy ampli- 
tude. In particular, we are going to discuss the solution of the equations 
for the asymptotic behaviour of the scattering amplitude at high energies. 

In the fifth section we will consider a general case of QCD and will dis- 
cuss the probabilistic interpretation as well as the scenario for the solution 
using the experience with the toy model. 

In conclusion we are going to compare our approach with other ap- 
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proaches on the market. 



2. The BFKL Pomeron Calculus 

The main ingredient of this calculus is the Green function of the BFKL 
Pomeron describing the propagation of a pair of gluons from rapidity y' 
and points x[ and x' 2 to rapidity y and points x\ and x 2 a . Since the 
Pomeron does not carry colour in the t-channel, we can treat initial and final 
coordinates as coordinates of quark and antiquark in a colourless dipolc. 
This Green function is well known 21 and has a form: 

G(x 1 ,x 2 ; y\x[,x' 2 ; y') = 0{y-y')x (1) 



OC 

X 

71— — OO 



^ / ^ i d2x o^ {n ' v){y ~ y ' ) Kn^)E(x u x 2 -,x \^E%x' 1 ,x' 2 -,x \u) 

= — oo 

where vertices E are given by 



E( Xl ,x 2 ;xoW) = \ -^ L -) ( -Phr J 5 ( 2 ) 

Xik = x i — Xk, Xi — Xi x + ixi y k ,x* — Xi x + ixi y ; h = (1 ^)/2 + iv and 
h = 1 — h*. The energy levels ui(n, v) are the BFKL eigenvalues 

w(n, v) = a s f - Re^-^P- + iv)] (3) 



where ^(z) = d\nT{z)/dz and Y{z) is the Euler gamma function, as = 
asNc/n and finally 

A( "' V) = [(n + l) 2 + 4^2]p - 1)2 + 4i ,2] ' ( 4 ) 

The interaction between Pomerons is described by the triple Pomeron 
vertex which can be written in the coordinate representation 5 for the fol- 
lowing process: two gluons with coordinates x\ and X3 at rapidity y\ decay 
into two gluons pairs with coordinates x" 1 and x" 2 at rapidity y 2 and x" 3 



a Coordinates Xi here are two dimensional vectors and, strictly speaking, should be de- 
noted by Xi or x;. However, we will use the notation Xi hoping that it will not cause 
difficulties in understanding. 

b Xi }X and Xi }V are components of the two dimensional vector Xi on a;-axis and y- axis 
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and x" 4 at rapidity y 3 due to the Pomeron splitting at rapidity y. It looks 
like 

T ,2 r Ai„i rfi ~_ rft n 



nai f cr x\ d 2 x 2 d 2 x 3 . , , , ,_. 

/ T 2 ^ r 2 {G(x 1 ,x 3 ;y 1 \x 1 ,x 3 ;y)^ 1 ^ 3 ,)- (5) 

Jv c J J '23 x 13 

G(xi , x 2 ; y I x" i , x" 2 ; y 2 ) G(4 , 4 ; y \ x" 3 , x" 4 ; y 2 ) 



where 



*&1',3' = r 'l3Pl>P3> Withp^ = -V 2 (6) 



and the arrow shows the direction of action of the operator L. For the 
inverse process of merging of two Pomerons into one we have 

7T a% f d 2 x[ d 2 x 2 d 2 x 3 „ , , „ „ , , 

^T" / 2 — 2 — 2 ^ i' 31 2;j/2|a; 1 ,a; 2 ;j/)G(a;"3,a; 4 ;y 2 \x 2 ,x 3 ;y) ■ 

JVc J X 12 X 23 X 13 

{^i',3'G{x'i,x' 3 ;y\xi,x 3 ;yi)) (7) 

The theory with the interaction given by Eq. (5) and Eq. (7) can be 
written through the functional integral 5 

= J D<f>D<P+e s with 5 = So + Si + S E (8) 

where So describes free Pomerons, Si corresponds to their mutual inter- 
action while Se relates to he interaction with the external sources (target 
and projectile). From Eq. (5) and Eq. (7) it is clear that 

s o = J dydy' d 2 x 1 d 2 x 2 d 2 x' 1 d 2 x' 2 (9) 
x 2 ; y) G~ 1 (xi,x 2 ; y\x[,x 2 ; y') ^(x[,x 2 ;y') 



irat 



5| f f d 2 xi d 2 x 2 d 2 x 3 



Sr = ^T I dV I J „2 '2 • (10) 



N c J J x 12 x 23 x 13 



■{(^ lt3 ^(x 1 ,x 3 ;y)) ■ $ + (x 1 ,x 2 ;y)<f> + (x 2 ,x 3 ) + h.c.} 

For S e we have local interaction both in rapidity and in coordinates, 
namely, 

S E = - J dyd 2 x Y d 2 x 2 (11) 
{$(x 1 ,x 2 ;y)T pr (x 1 ,x 2 ;y) + ® + (x 1 ,x 2 ;y)T tar (x 1 ,x 2 ;y) 
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where T pr (T tar stands for the projectile and target, respectively. The form 
of function r depends on the non-perturbative input in our problem and 
for the case of nucleus target they are written in Ref. 5 . 

For the case of the projectile being a dipole that scatters off a nucleus 
the scattering amplitude has the form 

where the extra as comes from our normalization. 

For further presentation we need some properties of the BFKL Green 
function 21 : 

1. Generally, 

G- 1 (x 1 ,x 2 ;y\x' 1 ,x' 2 ;y') = pjp 2 2 (^+#) = {j^ + H+ )p^ 

(13) 



as f d 2 xz x 2 
2ir J x 23 x 2 ^ 



12 



Hf(xi,xr,y) = (14) 



if(xi,x 2 ;y) - f{xi,x 3 ;y) - f(x 3 ,x 2 ;y)) 



2. The initial Green function (Go) is equal to 



. , , / - / 71-2 , X \ \iX\ y x \ v x\ 2 , 
Go{xi,x 2 ;y\x 1 ,x 2 ;y) = — In— — ln-^ — ^— (15) 

L X-y 2'X\i 2 2-1,2 X\i 2' 



3. It should be stressed that 

V 2 1 V 2 2 G (x u x 2 ;y\x' 1 ,x' 2 ;y') = (16) 

S (2) (x 1 -x' 1 )S (2 \x 2 -x' 2 ) + 5 (2) {x 1 -x' 2 )5 {2) {x 2 -x' l ) 

4. In the sum of Eq. (1) only the term with n = is essential for the 
high energy asymptotic behaviour since all u)(n, v) with n > 1 are negative 
and, therefore, lead to contributions that decrease with energy. Taking 
into account only the first term, one can see that G is the eigenfunction of 
operator L\ 3 , namely 

Li3G (x 1 ,X2;y\x' 1 ,x 2 ;y') = . Go{x 1 ,x 2 ;y\x' l ,x' 2 ;y') 

A(U, V) 

w G () (x 1 ,x 2 ;y\x / 1 ,x' 2 ;y'): (17) 
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The last equation holds only approximately in the region where v <C 1, but 
this is the most interesting region which is responsible for the high energy 
asymptotic behaviour of the scattering amplitude. 

Using Eq. (5), Eq. (7) and Eq. (2) we can easily obtain the chain equa- 
tion for multi-dipole amplitude T( n ) noticing that every dipole interacts 
only with one Pomeron (see Eq. (12)). The equation is shown in Fig. (1) 
and it has the following form 

3y = 2T J d K ^ x * z ) 

(T^(y;x 1 ,z) + T^(y;z,x 2 ) - T^\y;x 1 ,x 2 - T< 2 > (y; x u z; z, x 2 )) 

dT^(y;xi,x 2 ;x 3 ,x 4 ) a s f z , . , . 
dy = 2T J d K ^ x * z ) ^ 

(T^ ) (y;xi,z;x 3l x i ) + T (2) (y; z, x 2 ; x 3 , x 4 ) - T (2) (y; X\, x 2 ; x 3 , x 4 ) - 
-T^\y;xi,z;z, x 2 ;x 3 ,x 4 ) + {xi -> x 3 ,x 2 -> x 4 }) + 



+ ( ^T )2 ff / T ^i(^^x 3 ,x4x\xnV 2 x ,W 2 x ,T^(y-x',x''y, 

Deriving Eq. (18) and Eq. (19) we use Eq. (15) and Eq. (16) as well as the 
normalization condition (see Eq. (12)) for the scattering amplitude. These 
two equations are the same as in Ref. 22 . This shows that in the papers 
19 ' 22 , actually, the same approach is developed as in the BFKL Pomeron 
Calculus (much later, of course), in spite of the fact that the authors think 
that they are doing something more general. 

Assuming that = T {1 \ we obtain the Balitsky-Kovchegov 

equation 23 ' 24 . We can do this only if we can argue why the Pomeron 
splitting is more important than the Pomeron merging. For example this 
assumption is reasonable for the scattering of the dipole with the nucleus 
target. Generally speaking, the splitting and merging have the same order 
in as ( sec Eq. (5) and Eq. (7). In Eq. (18) and Eq. (19) these two processes 
look like having a different order of magnitude in as but this fact does not 
interrelate with any physics and reflect only our normalization. However, 
we will see that for a probabilistic interpretation the correct normalization 
is very important. 
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- 1 P -> 2 ^ 




+ ( 2P->P) -2 \P->2F\ : 



A3) 



Figure 1. The graphic form of equations for the multi dipole amplitude. 



Kernel K{x,y\z) is denned as 
K(x,y\z) = 
Function T^^i is equal to 



(x - y) 2 



(x — z) 2 (z — y) 2 
^2^\{x l ,y l - 1 x 2 ,y2\x,y) = 



(20) 
(21) 



d 2 zK(x, y; z) G (xi, yi\x, z) G (xi, yi\z, y) 



where Go is given by Eq. (15) and we use the notation (x, y) for a colourless 
pair of gluon (colour dipole). 

3. Generating functional and probabilistic interpretation 

In this section we discuss the main equations of the BFKL Pomeron Cal- 
culus in the formalism of the generating functional, which we consider as 
the most appropriate technique for the probabilistic interpretation of this 
approach top high energy scattering in QCD. 

To begin with, let us write down the definition of the generating func- 
tional 13 



Z(Y - Y ; M) 



(22) 



E 



/n 
Pn{Y - Y \ xi,yi;...;x i ,y i ]...i j x n! y n ) u{x l ,y i ) d 2 x^ 2 y t . 
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where u{xi,yi) = Ui is an arbitrary function of Xi and yi. The coordinates 
(xi,yi) describe the colorless pair of gluons or a dipole. P n is a probability 
density to find n dipoles with the size Xi—yi, and with the impact parameter 
(xi + yi)/2. It follows from the physical meaning of P n and the definition 
in Eq. (22) directly that the functional (Eq. (22)) obeys the condition 

Z(Y - Y ; [u = l]) = 1. (23) 

The physical meaning of this equation is that the sum over all probabilities 
is equal to unity. 

Introducing vertices for the dipole process: 1 — > 2 {V\^ 2 {x, y — ► X\,y\ + 
x 2 ,y 2 )), 2 —> 1 {V 2 ^i{xi,yi + x 2 ,y 2 —> x,y)) and 2^3 (V2— 3(^1) yi + 
x 2 ,y 2 —* x' 1 ,y' 1 +x 2 ,y 2 + x' 3 ,y' 3 ) we can write a typical birth-death equation 
in the form 

d P n (- ■ ■ ) Xj, yi; . . . ) _ . . 

dY - (M) 

= V 1 ^ 2 §Z){P n - 1 {...,x i ,y i ;...) - /•„; '•,.//,:...:•) (25) 

i 

X^-i (P n +i(;Xi,yi;...;x k ,y k ;) - P„(;x 4 ,y 4 ;)) (26) 

i>k 

^ V2^3^(Pn-i(;x i ,y i ;;x k ,yk;) ~ Pn(;xi,yi;;x k ,yk;)) (27) 



i>k 



Multiplying this equation by the product n™=i u i an d integrating over 
all Xi and yi, we obtain the following linear equation for the generating 
functional: 

dZ\u] 



dY 

with 



= X[u) Z[u] (28) 



X [u] = - J d i qd i q 1 d i q 2 (29) 

V 1 ^ 2 (q^q 1 +q 2 ) (-u(q) + u( qi )u(q 2 )) — — - (30) 

bu(q) 

-V 2 ^{q 1+ q 2 ^q) (u( qi ) u(q 2 ) u(q)) \ J^ffifoffi) " ( 31 ) 

Trying to make our presentation more transparent, we omitted in 
Eq. (29) the term that corresponding to the 2^3 transition (see Rcf. 
16 for full presentation). 
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Eq. (28) is a typical diffusion equation or Fokker-Planck equation 12 , 
with a diffusion coefficient depending on u. This is the master equation of 
our approach, and the goal of this talk is to find the correspondence of this 
equation with the BFKL Pomeron Calculus and the asymptotic solution 
to this equation. In spite of the fact that this is a functional equation 
we intuitively feel that this equation could be useful since we can develop 
a direct method for its solution and, on the other hand, there exist many 
studies of such an equation in the literature ( sec for example Ref. 12 ) as well 
as some physical realizations in statistical physics. The intimate relation 
between the Fokker-Planck equation and the high energy asymptotic was 
first pointed out by Weigert 25 in JIMWLK approach 27 , and has been 



= - T,n=i fln(xi,yi; ...;x n , y n ; Y ) n™ =1 -£-Z (Y, [m]) |„ i=1 d 2 x t d 2 y t 
= - Er=i( _1 )"/7n(a;i,yi; • ■ ■ ;x n ,y n ;Y Q )p(x 1 ,y 1 ;;x n ,y n ;Y - Y Q ) 



The physical meaning of functions j n is the imaginary part of the ampli- 
tude of the interaction of n dipoles with the target at low energies. All these 
functions should be taken from the non-pertubative QCD input. However, 
in Refs. 14 > 15 > 16 it was shown that we can introduce the amplitude of inter- 
action of n dipoles (7„(xi, yi; . . . ; x n , y n ; Y) at high energies (large values 
of rapidity Y) and Eq. (22), Eq. (28) and Eq. (32) can be rewritten as 
a chain set of equation for r ) n {x\, yi; . . . ; x n , y n ; Y). The equation has the 
form c 



discussed in Refs. 26 ' 19 > 20 . 

The scattering amplitude can be defined as a functional 24 ' 



N{Y;[li]) = 



(32) 



d-fn (ri,6i, . . . 



dY 





d A q'd i qV 1 ^ 2 (q l] qq') ln (...q' ...) 



c This equation is Eq. (2.19) of Ref. 16 but, hopefully, without misprints , part of which 
has been noticed in Ref. 22 . 
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n p n—1 p 

- £ / ^dYa^a^^™ (■■■,*■■■) -E MVV (33) 

n „ 

Vi-^fe; 9?')7n+i (■■■<!■ ••?')- E / rf4( ?^2^i g) 

• 7„_i (qt . . .qj . . .q) + d 4 qd 4 q'V 2 ^i (qqf, q') 7„_i (. . . ft . . . g) + 

i=i ■> 

n ,. 

+ / d4 Q v ^i(nqk;q) 7« (■••%• ••%•••) 

Comparing this equation for 71 = and 72 = T^ 2 ) one can see that 
V^ 2 = ^r^ 2 = ; (34) 

Z 7T Z 7T 

t/ /47r 2 a s 2 a s / r / , n , 

V2-1 = (— — ) — (-r 2 _ > i(a;i,j/i + Z2,y 2 z,y) + 

iV c Z 7T 

+ J d 2 xd 2 yT 2 ^ 1 (xi,yi+ x 2 ,y 2 ^ x,y)- 
■ (V 2 > ( Xl - X ) 6 W ( yi _ y ) + fiW (x 2 - x)5^ (y 2 - y))) (35) 
with r 2 _>i is given by Eq. (21). 



4. A toy model: Pomeron interaction and probabilistic 
interpretation 

In this section we consider the simple toy model in which the probabilities 
to find n-dipoles do not depend on the size of dipoles 13 ^ 14 . 16 . 17 . T n this 
model the master equation (28) has a simple form 

dZ dZ d 2 Z 

— = -r(l - 2) u(l ~u)^+ r(2 - 2) u(l -u)j^ (36) 

and this equation generates the Pomeron splitting Gp^ 2 p = —T(l — ► 2), 
Pomerons merging G 2 p^p = T(2 — > 1) and also the two Pomeron scattering 
G 2 p^ 2 p = — r(2 — > 1). It is easy to see that neglecting the u 2 d 2 Z/(du) 2 
term in Eq. (36), we cannot provide a correct sign for Pomerons merging 
G 2 p^p. 

Eq. (36) is the diffusion with the u dependence in the diffusion coeffi- 
cient. For u < 1 the diffusion coefficient is positive and the equation has 
a reasonable solution. If u > 1, the sign of this coefficient changes and 
the equation gives a solution which increases with Y and Z(Y) cannot be 
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treated as the generating function for the probabilities to find n dipoles 

(Pomerons) (see Refs. 9 ' n > 17 for details). We can see the same features in 

the asymptotic solution that is the solution to Eq. (36) with the l.h.s. equal 

to zero. It is easy to see that this solution has the form 

tvi ) 2 N 2 

Z{u;Y -oo) = e" with k = = > 1 (37) 

T(2 -> 1) a% 

One can see that for negative k this solution leads to Z > 1 for u < 1 . This 
shows that we cannot give a probabilistic interpretation for such a solution. 

Using the asymptotic solution and finding the typical u at high energies 
from the unitarity constraints 17 > 28 > 29 we can determine the asymptotic limit 
for the high energy amplitude. It turns out that 

N(Y -> oo) = 1 - e~ K (38) 

It was shown in Ref. 17 that we can search for the correction to the 
solution of Eq. (38) in the form 

Z(u;Y) = Z(u;Y — oo) + AZ(u;Y) (39) 

Asssuming AZ(u; Y)(u; Y — ► oo) ^> Z(u; Y — > oo) we can write for 
AZ(u; Y) the linear equation 17 and the solution of this equation decreases 
with energy 

It is interesting to notice that this asymptotic amplitude does not show 
a black disc behaviour at high energies. The behaviour of the scattering 
amplitude, given by Eq. (38), corresponds to the so called gray disc be- 
haviour. 

5. Probabilistic interpretation in QCD 

Eq. (34) has a very simple physical meaning describing the Pomeron split- 
ting as the decay process of one dipole into two dipoles. It turns out that 
the vertices for 2P — ► 3P and 2P — ► 2P can be easily understood as a 
dipole 'swing'. What we mean is that with some probability two quarks 
of a pair of dipoles can exchange their antiquarks to form another pair of 
dipoles 16 . Naturally, this process has an 1/N 2 suppression and it correctly 
reproduces the splitting and rescattering of two Pomerons that has been 
explicitly calculated from the diagrams 7 . 

Eq. (35) d is more difficult to view as the vertex for the transition of 
two dipoles into one. Indeed, the integral over coordinates of the produced 

d This equation is quite different from the equation which is obtained in Ref. 16 (see also 
Ref. 22 ). The main difference stems from the correct use of the BFKL Pomeron Calculus 
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dipole is positive, namely / V(2 — » l)d 2 xd 2 y = J T(2 — > l)d 2 xd 2 y. How- 
ever, Eq. (35) leads, generally speaking, to a negative vertex in some regions 
of the phase space. The detailed study of this vertex will be published soon 
18 . Here we want to point out that the key problem is not in the probabilis- 
tic interpretation of the microscopic process of two dipoles to one dipole 
transition but the fact that a negative vertex V(2 — > 1) means that in some 
kinematic region we have a negative diffusion coefficient which results in a 
solution that increases at large values of rapidity Y. To save such a theory, 
we have to introduce other Pomcron interactions like 2P — > 3P and/or 
2P — > 2P transitions. 

In Ref. 16,17 attempts were made to deal with such Pomeron interac- 
tions. It turns out that at large values of Y the 2 — > 1 process contribute 
in a very limited part of the kinematic region with a very specific function 
u(x, y) . In this particular region the vertex V{2 — ► 1) is positive 18 . There- 
fore, we could use the probabilistic interpretation but we need to study this 
process better and deeper to obtain the final result. 

In the toy model it has been shown 9,11 that we can generate the 2P — > 
3P vertex without the process of two dipoles to one dipole transition. The 
transition of two dipoles to two or more dipoles also leads to this vertex. 
Similar ideas are developed in Ref. 31 in QCD. However, we need to pay 
a price: the contribution to Eq. (28) will be negative to give the correct 
sign for 2P — > P interaction. In other words, we can add to Eq. (28) the 
contribution 

.2 2 

/ W d 2 Xi d 2 y t JJ d 2 x'i d 2 y[ V 2 ^ 2 {x x ,yi + x 2 ,y 2 -> x[, y[ + x' 2 ,y' 2 ) 
J i=l i=l 

i , /n / / /x 1 d d 

u{x 1 ,y 1 )u(x 2 ,y 2 ) -— — 

2du(x 1 ,y 1 ) du(x 2 ,y 2 ) 

which will give the 2P —> P vertex in the form 

T(2P^P)(x 1 ,y 1 +x 2 ,y 2 ^x' 1 ,y' 1 ) = 

- J d 2 x' 2 d 2 y' 2 V 2 ^ 2 {x x ,y x + x 2 ,y 2 x\,y\ + x' 2 ,y 2 ) 

Therefore, we have to assume that V 2 ^ 2 < 0. As we have discussed, 
generally speaking, it means that our problem has no solution. However, in 



for determining this vertex while in Refs. 16 > 22 the Born diagram was used which docs 
not and cannot give a correct expression. 



(40) 
(41) 
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QCD the situation is much better. Indeed, to generate a correct 2P — > P 
vertex we need to introduce V2— >2 °c The Feymann diagrams for this 
transition are obvious (see Fig. (2-b). However, the main contribution 
stems from the diagrams of Fig. (2-a) - type which are of the order of 
a$/N^. Therefore, diagrams of Fig. (2-b) - type are small corrections to 
the main contribution and could be negative, in spite of the fact that the 
diagrams shown in Fig. (2-b) actually give a positive contribution. We will 
discuss this scenario in all details in our paper 18 . 




6. Conclusions 

Being elegant and beautiful the BFKL Pomeron Calculus has a clear disad- 
vantage: it is not a theoretically closed theory. Indeed, we need to add to 
the formalism of the BFKL Pomeron Calculus the theoretical ideas what 
kind of Pomeron interactions we should take into account and why. Of 
course, the Feymann diagrams in leading log l/x approximation of per- 
turbative QCD allow us, in principle, to calculate all possible Pomeron 
interactions but, practically, it is a very hard job. Even if we will calcu- 
late these vertices we need to understand what set of vertices we should 
take into account for the calculation of the scattering amplitude. This is 
the reason why we need to develop a more general formalism. Fortunately, 
such a formalism has been built and it is known under the abbreviation 
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JIMWLK-Balitsky approach 30 > 27 > 23 . I n this approach we are able to cal- 
culate all vertices for Pomeron interactions as it was demonstrated in Ref. 
31 and it solves the first part of the problem: determination of all possible 
Pomeron interactions. However, we need to understand what vertices we 
should take into account for the calculation of the scattering amplitude. We 
hope that a further progress in going beyond of the BFKL Pomeron Calcu- 
lus (see Refs. 31 < 32 ) will lead to such a development of the BFKL Pomeron 
Calculus that we will have a consistent theoretical approach. Hopefully 
this approach will be simpler than Lipatov's effective action 33 which is not 
easier to solve than the full QCD Lagrangian. 

The probabilistic interpretation gives a practical method for creating a 
Monte Carlo code in spirit of the approach suggested in Ref. 34 . This code 
will allow us to find a numerical solution to the problem and to consider 
the inclusive observables. This extension is very desirable since most of the 
experimental data exist just for these observables. 

This year my teacher, Prof. Gribov, would have been 75 years old. I am 
happy to report him that his ideas are still alive and still give the simplest 
and the most elegant approach to high energy interaction. The high energy 
behaviour of the scattering amplitude is still an open theoretical question. 
We have learned a lot and we are able to foresee a lot of difficulties ahead. 
Therefore, it is a good problem to be solved. 
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